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The generalized time-dependent harmonic oscillator is studied. Though several 
approaches to the solution of this model have been available, yet a new approach is 
presented here, which is very suitable for the study of cyclic solutions and geometric 
phases. In this approach, finding the time evolution operator for the Schrodinger 
equation is reduced to solving an ordinary differential equation for a c-number vector 
which moves on a hyperboloid in a three-dimensional space. Cyclic solutions do not 
exist for all time intervals. A necessary and sufficient condition for the existence 
of cyclic solutions is given. There may exist some particular time interval in which 
all solutions with definite parity, or even all solutions, are cyclic. Criterions for 
the appearance of such cases are given. The known relation that the nonadiabatic 
geometric phase for a cyclic solution is proportional to the classical Hannay angle 
is reestablished. However, this is valid only for special cyclic solutions. For more 
general ones, the nonadiabatic geometric phase may contain an extra term. Several 
cases with relatively simple Hamiltonians are solved and discussed in detail. Cyclic 
solutions exist in most cases. The pattern of the motion, say, finite or infinite, can not 
be simply determined by the nature of the Hamiltonian (elliptic or hyperbolic, etc.). 
For a Hamiltonian with a definite nature, the motion can changes from one pattern 
to another, that is, some kind of phase transition may occur, if some parameter in 
the Hamiltonian goes through some critical value. 

PACS numbers: 03.65. Ca, 03.65.Vf 

I. INTRODUCTION 

The harmonic oscillator is one of the most familiar models in physics. It is widely used in 
various fields, both classically and quantum-mechanically. Several physical problems can be 
described by the harmonic oscillator with time-dependent parameters [1, 2]. More generally, 
such time-dependent parameters may describe approximately the interaction of the harmonic 
oscillator with some external degrees of freedom. Therefore, the time-dependent harmonic 
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oscillator has been the subject of many theoretical researches over decades [3-6]. Since the 
discovery of the geometric phase [7-14], the model has attracted more attention [15-22], be- 
cause it is simple and serves as a good example for the study of geometric phases, just like the 
case of particles with spin and magnetic moment moving in time-dependent magnetic fields 
[23-37]. The two models are similar in that they involve similar Lie algebras. In fact, the 
Hamiltonian for spin in a magnetic field is an element of the S0(3) algebra while that for 
the time-dependent harmonic oscillator is one of the S0(2,l) algebra. However, the time- 
dependent harmonic oscillator is of more interest since it has a classical counterpart, and the 
quantum motion can be compared with the classical one. For example, the relation between 
the geometric phase and the classical Hannay angle [38-40] is an interesting subject. 

There exist mainly two approaches to the solution of the time-dependent harmonic oscil- 
lator. The main point of the first approach [3, 4] is to find an invariant operator and its 
eigenstates. The second approach employs time- dependent unitary transformations [6, 16]. 
Here we develop another approach to the solution of the Schrodinger equation. It is a fur- 
ther development of the approach previously used for spin moving in time-dependent magnetic 
fields [36]. First we find an invariant operator, and then go further to obtain the time evolution 
operator. It should be remarked that our method of finding the invariant operator is rather 
different from that in the first approach. In that approach the problem is reduced to solving 
a nonlinear differential equation. In our approach, it is reduced to solving a linear differential 
equation for a three- component c-number vector, and thus is simpler. On the other hand, the 
time-dependent unitary transformation approach in Ref. [16] seems even simpler than ours, 
but our approach can be easily generalized to other systems where the Hamiltonian is an ele- 
ment of a more complicated Lie algebra. However, the main advantage of our approach is that 
it is very suitable for the study of cyclic solutions and geometric phases. 

In the literature there exists some argument that cyclic solutions are available for any time 
interval, say, [0, r] where r is arbitrary, because one can always choose the eigenstates of U{t), 
the time evolution operator at r, as initial conditions at t = 0. This is true. However, the 
problem is that, for the time-dependent harmonic oscillator, U (r) may have no normalizable 
eigenstate. This is different from the case of spin, where no problem of normalizability has to be 
worried. Obviously, normalizable states are physically more interesting than nonnormalizable 
ones. Moreover, if one is interested in the geometric phases, it seems still not clear how to 
define it for the nonnormalizable cyclic solutions. The difficulty lies in the definition of the 
dynamical phase. If only normalizable states are to be considered, then one should be able to 
tell whether there are cyclic solutions for a given r. This is not a trivial task even when the 
time evolution operator is explicitly available. In our approach, however, the problem can be 
solved naturally. We will give a necessary and sufficient condition for the existence of cychc 
solutions in an arbitrarily given time interval. 

It has been shown by several authors that the nonadiabatic geometric phase for a quantum 
cyclic state is equal to —{n+ 1/2) times the classical Hannay angle [19-21]. However, it 
seems not very clear under what restriction on the initial condition is this relation valid, or 
whether modification to this relation is needed in some case where cyclic solutions exist for 
less restricted initial conditions. We will reestablish the above relation and show that it is 
valid only for cyclic solutions with special initial conditions. There exist several cases where 
more general cyclic solutions exist in some particular time interval. Among these cases two 
are of special interest. In one of the two cases all solutions are cyclic, and in the other case all 
solutions with definite parity are cyclic. We will give criteria for the appearance of such cases. 
In all these cases, the nonadiabatic geometric phase contains in general an extra term which 
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depends on the initial condition, in addition to the one proportional to the classical Hannay 
angle. Similar situations have been found in other systems [35, 36, 41]. 

A normalizable state can be regarded as a wave packet in the configuration space. The 
motion of a wave packet can be roughly described by the change of its position and width. For 
the time-dependent harmonic oscillator, we will show that if the position of a wave packet is 
confined in a finite region, then its width also remains finite, and the contrary is also true. A 
problem that seems not clear concerns the relation between the pattern of the motion, say, finite 
or infinite, and the nature of the Hamiltonian, that is, elliptic, hyperbolic, or critical (see Sec. 
II). In simple cases where the Hamiltonian is time independent (or the time dependence lies 
only in an overall factor) , elliptic Hamiltonian leads to finite motion and the other ones leads to 
infinite motion. However, if the Hamiltonian is time dependent, the situation is complicated. 
We will see that for a Hamiltonian with a definite nature, say, elliptic, the motion of the wave 
packet may exhibit different patterns if some parameter in the Hamiltonian takes different 
values. In particular, when the parameter goes through some critical value, some kind of 
phase transition occurs, that is, the motion changes from one pattern to another. At the 
critical value, the motion has an independent pattern. 

This paper is organized as follows. In Sec. II we develop some mathematical formulas that 
will be used in the subsequent sections. In Sec. HI a new method is presented to derive the 
time evolution operator for the Schrodinger equation. In Sec. IV a necessary and sufficient 
condition for the existence of cyclic solutions in an arbitrarily given time interval [0, r] is 
given, and the known relation between the nonadiabatic geometric phase and the classical 
Hannay angle is reestablished. In Sec. V we study several cases where more cyclic solutions 
are available, and give criterion for the appearance of such cases. A modification to the above 
relation between the nonadiabatic geometric phase and the classical Hannay angle are also 
discussed in this section. In Sec. VI we study several examples where explicit results are 
available. The existence of cyclic solutions are discussed in detail, and it is seen that they exist 
in most of the examples. The evolution of normalizable states are also studied from the point 
of view of wave packets. Various patterns of motion are revealed in these examples, and phase 
transition is explicitly observed. Sec. VII is devoted to a brief summary. In the appendix we 
briefiy discuss how to extend our formalism to a more general system. 



II. THE MODEL AND SOME MATHEMATICS 

The time-dependent harmonic oscillator is described by the Schrodinger equation 

ihdtilj{t) = H{t)ilj{t) (la) 

with the following Hamiltonian. 

H{t) ^u{t)K -n^it), (lb) 

where uj{t) is a time- dependent frequency parameter, n(t) = {ni{t) , n2{t) , n^{t)) is a time- 
dependent c-number vector and K = {Ki, K2, K^) is an operator vector defined below. In 
this paper, any vector a — (01,02,03) has an associated vector — (—01,-02,03), and 
scalar product between two vectors a and b always appears as (a, 6) — a • — ■ b — 
O363 — O161 — 0262- In matrix form these are = ga and (a, b) = a^gb = b^ga where a and b are 
now column vectors and the superscript t denote transposition, and g = diag(— 1, —1, 1). In 
the following the square of a vector always stands for a ■ rather than a - a. In the above 
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Hamiltonian we always choose uj such that equals 1, —1 or 0. Thus among the components 
of n only two are independent, and the Hamiltonian involves three independent parameters. 
The vector K is defined as 

i^i - J f /xoa;o^' - , K^^-UxP + PX), ^3 = i ("^xoc^oX^ + ^ y (2a) 

2 V A*o<^o/ 2 2 V A*o<^o/ 

where X is the coordinate and P the momentum, satisfying [X, P] — ih, /iq is the mass of the 
particle, and uq is some constant frequency parameter. For simplicity in notations, we define 
X = ^fiQiVoX and p = P/ y/fiQUjQ, which still satisfy the commutation relation [x, p] — ih, then 
the above expressions become 

Xi = |(x2-p2), K2^-l(xp + px), Ks^lix^+p-"). (2b) 

Note that both x and p have the same dimensionality as ^/h, and K has the same dimension- 
ality as h. The components of K satisfies the commutation relation 

[Ki,Kj]^-2iheijkKl. (3) 

An equivalent form is [Kf, K^] = —2iheijkKk. This is an S0(2,l) algebra. The Hamiltonian is 
said to be elliptic, hyperbohc or critical if equals 1, —1 or 0, respectively. 

To solve the Schrodinger equation we need some operator formulas. We will briefly derive 
them here. Let 

F{0 = b)KQ^{^, b), g(e, b) = exp(-iin-ie^ ■ b^), (4) 

where ,^ is a real parameter and b a real vector independent of ^. This is a unitary transfor- 
mation of K. Using Eq. (3) it can be shown that 

F'(0 = 6^xF^(e), (5) 

where the prime indicates derivative with respect to ^. An equivalent equation is F^'{^) = 
b X F{^), which is convenient in obtaining 

F"{0 + b'F{0^[F{0-b^]b. (6) 

However, Eq. (5) leads to F'(C) ■ 6^ = 0, and thus F(0 • &^ = F{0) ^ K l^. Therefore 
Eq. (6) is simplifled as 

F"(0 + 6^F(0 = {K . ¥)b. (7a) 
This is a simple equation. With the initial condition 

F(0) = K, F'(0) = ¥ X K^, (7b) 

the solution is easily flxed. The results are listed below. 
If 6^ = 1, we denote 6 by 6+, and have 

g(C, b+)KQ^{i, 6+) = [K-{K- bl)b+] cose + bl x K^sin^ + {K ■ 6^)6+. (8) 

When e = 2iV7r where is an integer, we have Q{2N7r,b+)KQ^{2N7r,b+) = K, or 
Q{2N7r,b+)K = KQ(2NTi,b+). That is, Q(2NTi,b+), and in particular, exp(i;i~iA^7rA'3), 
commutes with K. However, this does not mean that Q{2Nt:, 6_|_) is a c-number, since a 
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c-number must commute with x and p. But we will see below that Q{AN'k, 6+) is indeed a 
c-number. 

If 6^ = —1, we denote 6 by 6_, and have 

Q{C,b-)KQ\C,b-) = [K + {K ■yL)b-] cosh ^ + b^ x sinh^ - (K • 6^)6_. (9) 

When 6- = (— sin 0, cos0, 0), the result is useful in the following sections. In this case we 
denote 6_ by 6^ and Q{^,b-) by 0): 

Q(e,(/.) = exp{-lih-^CK ■ bl) = exp[-iin-ie(i^isin0 - Kscos^)]. (10) 

In particular we write down 

g(^,0)/s:3Qt(^,0) = /sTg cosh ^ - /sTisinh^ cos 0-i^2sinh^ sin 0. (11) 

If 6^ = 0, we denote 6 by feo, and have 

Q(e, bo)KQ\^, b,) = K + xK^ + ie{K ■ 6g)6o. (12) 

For a c-number vector a, one can define a transformed vector A by 

K-A^Q{^,b)K-aQ\tb). (13) 

It can be shown by straightforward calculations that A"^ — a^, regardless of the values of ^ 
and b. Thus the elements of the S0(2,l) algebra, K • a, are distinguished into three classes, 
characterized by whether is positive, negative or zero. They cannot be connected by any 
one of the above unitary transformations. For example, K?, cannot be transformed to Ki and 
vice versa. Geometrically, the surface defined by = is a cone in the a space. Therefore, 
the three classes of elements are characterized by whether a is inside, outside or on the cone. 

Similarly, one can calculate the unitary transformation of x and p. Here we only write down 
the result for the case with — 1 (the subscript + in the components of 6+ is omitted) . 

(5(C, b+)xQ\^, b+) = X ^cos I - 62 sin - p{b^ + 61) sin |, (14a) 



g(e,6+)pgne,&+) =^(&3-?^i)sin^+p(^cos^ + 62sin^J . (14b) 

Here we see that Q{AN7r, 6+) commutes with x and p. Thus it must be a c-number, as 

mentioned above. 

The final point of this section concerns the eigenvalues and eigenstates of the operator K. 
We are only interested in normalizable, or bound states. Ki and K2 do not have normalizable 

eigenstates. The eigenvalues of K3 are knh where A;„ = n + 1/2 and n — 0,1,2, The 

corresponding eigenstates will be denoted by ipn- In the coordinate representation, ipn{x) — 
Nnexp{—x'^/2h)Hn{x/Vh), where are Hermite polynomials and = (l/2"n!-\/7r^)^/^. If 
a vector e satisfies > 0, then K e^ has normalizable eigenstates. Without loss of generality, 
we take — 1 and 63 > 0, then e can be written as 

e — (sinh^cos0, sinh^ sin 0, cosh ^). (15) 

According to Eq. (11), K ■ — Q{^, 4>)KsQ^{^, (f>). Therefore the eigenvalues oi K • are 
still knh, and the corresponding eigenstates are 

C = 0)^n = exp{-lih-'^K ■ b^^n. (16) 
This result will be employed below. 
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III. TIME EVOLUTION OPERATOR 

In this section wc deal with the time evolution operator for the Schrodinger equation (1). 
We define a time-dependent c- number vector e(t) by the differential equation 

e(t) = -2uj{t)n%t) X e%t) (17a) 

where the overdot denotes differentiation with respect to and the initial condition 

e(0) = eo, (17b) 

where Cq satisfies 

el = 1, eo3 > 0, (17c) 

and is otherwise arbitrary. We would assume that n{t) and uj{t) varies continuously, so that 
any solution e{t) is well behaved. If one solution to this equation can be found, then the time 
evolution operator for the Schrodinger equation can be worked out. It should be remarked 

that the above equation for e{t) is the one satisfied by {K), the mean value of the operator 
vector K in an arbitrary state (cf. Eq. (27) below). More discussions can be found in the 
appendix. 

We take the initial state of the system to be '0(0) = V'n"' ^^^^ 

K ■ elip{Q) ^ knhip{Q), n = 0,1,2,.... (18) 
If evolves according to Eq. (1) and e evolves according to Eq. (17), then 

K ■ e\t)il]{t) = knhiljit) (19) 

would hold at all later times. This can be easily proven by induction. 

By definition, Eq. (19) is valid at t = 0. We assume that it is valid at time t, what we need 
to do is to show that it is also true at time t -\- At where At is an infinitesimal increment of 
time. In fact, using Eqs. (1) and (17) we have 

V'(t + At) = V'(i) - ih-^ou{t)K ■ n%t)iP{t)At, (20a) 

e(t + At) = e(t) - 2uj{t)n\t) x e^{t)At. (20b) 

After some simple algebra, the conclusion is achieved. 

It should be remarked here that K ■ e^{t) is an invariant operator, so that it has time- 
independent eigenvalues. Indeed, it is easy to show that 

ihK ■ e^{t) + [K ■ e^{t), H] = 0. (21) 

We see that once a solution e{t) is found, an invariant operator is obtained. Since e{t) satisfies 
a linear differential equation, this method is convenient. Moreover, it can be easily generalized 
to systems where the Hamiltonian is an element of a more complicated Lie algebra. A brief 
discussion on this point is given in the appendix. 

From Eq. (17) it is easy to show that e^(t) = Bq = 1. This yields |e3(t)| > 1. As is assumed, 
e{t) varies continuously, thus 63 (t) keeps its original sign at all later times. Therefore, e{t) can 
be written in the form of Eq. (15), where ^ — C{t) and = 0(t), and 



iP{t) = exp[ian{t)]Q{i{t),<P{tM 



(22) 
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where Q;„(t) is a phase that cannot be determined by the eigenvahie equation. However, an{t) 
is not arbitrary. To satisfy the Schrodinger equation, it should be determined by the other 
variables ^{t) and 4>{t). In fact, the above equation yields 

(V'W, + At)) = 1 + idn{t)At + (V'n, QHt <t>)dtQ{^, (t>)ll^n)^t. (23) 

Using the formula [42] 

ent)d^^-F{t) ^ _ ^XF(f)p^^y-XFif) (24) 

where F{t) is any operator depending on then using Eq. (9), and notice that (■^n, i^iV'n) — 
(ijjn, K2il^n) = 0, we obtain 

{il){t),%l){t + /^t)) = l + id„(t)At- |ifc„0(t)[cosh{(t) - l]Ai. (25) 

On the other hand, from Eq. (20) we have 

iip{t),ij;{t + At)) = 1 - ih-^uj{t)u{t) ■ n'^{t)At, (26) 

where 

u{t) = (m,Km)- (27) 

This definition will be repeatedly used below. It is easy to show that u{t) satisfies the same 
equation as e(t), and for the above initial state u(0) = /c„/ieo, so we have u{t) — knhe{t). 
Comparing the two results above and taking this relation into account, we obtain 



Therefore 
where 



oin{t) = \kJ{t)[cos\ii{t) - 1] - Kuj{t)e{t) ■ n3{t). (28) 
an(t)-an{0)^kna{t), (29) 



a(t) = - /" <j){t')[cosh^(t') - 1] dt' - I uj{t')e{t') ■ n%t') dt' . (30) 
Substituting into Eq. (22) we obtain 

m = Qim, m) e^m-'a{t)K^]Q\m, m)m- (31) 

We denote the time evolution operator as U{t), defined by the equation il>{t) = U{t)ijj{0) with 
an arbitrary ■0(0), then the above equation is equivalent to 

um:° = Qim, ^w) exp[ih-Mt)Ks]Q^{m, m)^n°- (32) 

Now an arbitrary normalizable initial state ■0(0) can be expanded as 

^(0) = ^c„C°- (33) 

n 

Applying U (t) to both sides of this equation, using Eq. (32), and noting that the operators on 
the right-hand side of that equation is independent of n, we immediately realize that Eq. (31) 
is in fact valid for an arbitrary initial state. Thus we arrive at the result 

U{t) = exp[-lm-^^{t)K ■ blit)] exp[m-^a{t)K3] exp[hh-^^{0)K ■ 6^(0)]. (34a) 
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Using Eq. (9), it can be recast in the form 

U{t) = eM-l^^-'mK ■ eM^fi-'mK ■ 6^(0)] exp[in-ia(t)K • eg]. (34b) 

Eq. (34b) is suitable for the general discussions below while Eq. (34a) may be more convenient 
for practical calculations. 

Let us make some remarks on the result. First, we see that once a solution of Eq. (17) is 
found, the time evolution operator for Eq. (1) is available. The result depends formally on Sq, 
but Bq is merely an auxiliary object, hence the result must be essentially independent of it, 
though it might be difficult to prove this explicitly. On the other hand, it is the flexibility in the 
choice of Bq that makes it convenient for the general discussions of cyclic solutions. In practical 
calculations, one should choose a solution e{t) that is as simple as possible such that U{t) can 
be easily reduced to the simplest form. Second, the operator U{t) depends not only on e(t), 
but also on the history of it. This is obvious from Eq. (30). Third, though 0(t) is indefinite 
when C,{t) = 0, it is obvious that U{t) is well behaved everywhere. Fourth, by straightforward 
calculations it can be shown that ihdtU{t) = oo{t)K ■ n^{t)U{t) and f/(0) = 1, as expected. In 
other words, though U (t) is obtained by considering the evolution of normalizable states, it is 
also valid for nonnormalizable ones. 



IV. CYCLIC SOLUTIONS AND GEOMETRIC PHASES 

Now we can go further to discuss cyclic solutions in any time interval [0, r] where r is an 
arbitrarily given time. These cyclic solutions are not necessarily cyclic in subsequent time 
intervals with the same length, say, [r, 2r] . 

Since Eq. (17) is a linear differential equation, the general solution e{t) must depend on 
the initial vector Cq linearly. Thus it can be written in a matrix form 

e{t) = E{t)eQ, (35) 

where e{t) and eo are column vectors, and E{t) is a 3 x 3 matrix which is obviously real. If 
both e^'^\t) and e^'^\t) are solutions to Eq. (17), it is easy to show that e^^\t) ■ e^^^^{t) = 
e(i)(0) • e(2)s(0). Therefore the matrix E{t) satisfies 

E\t)gE{t)^g. (36) 

This yields det E{t) = ±1. As is assumed, E{t) varies continuously, and detE{0) = 1, so that 
det E{t) — 1. Therefore the product of the three eigenvalues of E(t) must be 1, and none can 
be zero. Now if e„ is an eigenvector with eigenvalue a, that is, Ee„ = ae„, it can be easily 
shown that E^{gecr) = (7~^{gecr). This means that a^^ is an eigenvalue of E^, and thus an 
eigenvalue of E. Therefore the eigenvalues of E{t) should be {l,a(t),a^^{t)}. Since E is real, 
a* is its eigenvalue if a is one. Therefore, if a{t) is complex, it must be unit: |cr(t)| = 1. 

If cr(r) 7^ 1 at the time r, one eigenvector 17 (r) of the matrix E{t) with eigenvalue 1 can 
be found, which satisfies E{T)r]{r) — r]{r). It can be taken as real. t]'^{t) may be positive, 
negative or zero, depending on E{t) and r. First we consider the case with 

V'{r) > 0. (37) 

Because 77 (r) is only determined up to a constant factor, we can choose that constant such 
that r;^(r) — 1 and 773 > 0. Then we can take 



eo = t7(t) 



(38) 
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as the initial condition in Eq. (17), and have e(T) = E{T)eo — E{T)r]{T) — r]{T) — Cq, that is 

e(r) = eo. (39) 

This means that ^(r) = ^(0) and ^^(t) = 6^(0), and leads to 

U{t) = ex:p[ifr^a{T)K ■ e^] . (40) 

Now it is clear that with the initial condition '(/'(O) = ip^ (n = 0, 1,2, . . .), we have a cyclic 
solution in the time interval [0,r]. More specifically, tpir) = e^^"ip{0), where the total phase 
change is 5„ = A;„Q;(T),mod 2%, with q;(t) given by 

(r) = i I 0(i)[cosh^(i) -l]dt- f uj{t)e{t) ■ n%t) dt. (41) 

For the present state, u{t) — knhe{t), so the dynamic phase /?„ = —h~^ {H{t)) dt turns out 
to be 

Pn = -kn [ cu{t)e{t) ■ n%t) dt. (42) 
Jo 

Therefore the nonadiabatic geometric phase jn — — Pn is given by 

7n = —kn^Og, mod 2n, (43) 

where 

.e,---^m^..m-.u.---^^^^ ,44) 

will be shown to be the classical Hannay angle below. Because e^(i) — 1 and 63 (t) > 0, e{t) 
moves on the upper sheet of a hyperboloid. This is a basic consequence of the fact that the 
Hamiltonian is an element of the S0(2,l) algebra. The above integral can be recast in two 
other forms 

A0 ^ - f d5 = - f d5i2 
' ^2Js^ei + ei + ei ^ 2 J s,, ^lT4 + 4^ 
where S is the surface enclosed by the closed trace of e{t) on the hyperboloid, and d^* the 
surface element; S12 is the projection of S on the 6162 plane, and dS'12 the area element; the 
upper (lower) sign corresponds to an anticlockwise (clockwise) trace of e{t). The geometric 
nature of the Hannay angle is obvious from the above expression, because it depends only on 
the closed trace of e(t), but not on the details of the traversing process. Because of the relation 
(43), the geometric nature of the nonadiabatic geometric phase is also obvious. 

Thus Eq. (37) is a sufficient condition for the existence of cyclic solutions. Under this 
condition there exist at least a denumerable set of normalizable cyclic solutions in the time 
interval [0,t]. Of course they may be trivial ones in some cases. All phases can be expressed 
in terms of the vector e{t). The relation between the nonadiabatic geometric phase and the 
Hannay angle is reestablished. 

States with initial condition other than the above ones are in general not cyclic ones, 
even though in these initial states u(0) points in the direction of Cq such that u{t) = u(0). 
However, if Q;(r)/7r happens to be a rational number, more cyclic solutions are available, and 
the above relation between the nonadiabatic geometric phase and the Hannay angle would 
need modification for these cyclic solutions. This will be discussed in the next section. 

If r7^(r) < 0, one can still take e(0) = t7(t) as the initial condition for Eq. (17a), and have 
e(r) = e(0). However, this solution cannot be used in the time evolution operator (34), and no 
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similar discussions to the above are available. In fact, there is no normalizable cyclic solution 
in this case, since the condition (37) is also a necessary one. This is proven below. 

If there exist one normalizable cyclic solution in the time interval [0,t], that is, iP{t) = 
e^^ip{0), then in this state u{t) = u{0), or u{t) = u{0) in the form of column vectors. As 
pointed out before, u{t) satisfies the same equation as e(t), thus u{t) = E{t)u{0). Comparing 
the two relations we obtain E(t)u{0) — u{0). In other words, u{0) is an eigenvector of E{t) 
with cigcnvahic 1. The remaining point is to show that w^(0) > 0. Indeed, for any normalizable 
state ip{t), it is not difficult to show that u'^it) > > 0, by using the Schwarz inequality. 

To conclude this section let us work out the Hannay angle in terms of e(t). We denote the 
classical coordinate by qc and momentum by p^, and define a vector / as 

h^W!-pl): h^-qcPc: h^\{ql+pl). (46) 

The classical Hamiltonian is now — oj{t)I ■n^{t), and the evolution of qc and Pc are governed 
by the canonical equations of motion. It is easy to show that I — I ■ e^{t) is an invariant. This 
leads to a quadratic equation in qc and Pc'- 

(es + ei)pl + 2e2gcPc + (eg - ei)g,' = 21. (47) 

If = 1, this describes an ellipse on the qcPc plane, whose area is 27r/. The ellipse changes 
its shape when e(t) varies with time. If e(r) = e(0), the ellipse at the time r coincide with 
that at the initial time. This is a classical nonadiabatic cyclic evolution. The qc and Pc can be 
expressed in terms of / and its canonical variables 9 as 



qc(9,I,e) = ^/2I(e3 + ei) cos 6', Pc(9,I,e) = -\ — — — (e2 cos6' + sin6'). 

V 63 + ei 



(48) 



Using = 1, we have dpcAdgc = /cos^ 6'dei Ade2/e3, and the contour average is (dpcAdgc) 
Idei A 6.62/263. According to Ref. [40], A^g = —dj f {dpc A d^c), we arrive at 



A9, 



_1 r de,Ade2_ 
2 Jsi2 ea 



Note that dci Ade2 corresponds to dS'12 (— dS'12) for an anticlockwise (clockwise) trace of e{t), 
and 63 = ^/l + el + el, this is the same as Eq. (45). It is independent of /. 



V. MORE ON CYCLIC SOLUTIONS AND GEOMETRIC PHASES 

In the last section we have shown that Eq. (37) is a sufficient and necessary condition for 
the existence of cyclic solutions in the time interval [0, r]. Under this condition, there exist 
a denumerable set of cyclic solutions. In this section we discuss some special cases where 
more general cychc solutions are available. We will see that the simple relation (43) has to be 
modified. 

Before discussing these cases, we define 

for any state ip{t), and study their evolution with time. According to the Schrodinger equation, 
they satisfy the equation of motion: 

X = uj[n2X + {ni + n3)p], p = uj[{ni - n3)x - n2p]. (51) 
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This is the same as that for the classical variables and p^-, since the Hamiltonian is quadratic 
in X and p. We denote a two-component column vector q — Because the about equation 

is linear, we have 

q{t) = E^{t)q{0), (52) 
where Eg{t) is a 2 x 2 evolution matrix independent of q{0). Next we define a vector v{t) as 

vi = l{x^ -p^), V2 = -xp, V3 = l{x^+p'^). (53) 

It is straightforward to show that v(t) satisfies the same equation of motion as e(t) or u{t). 
Therefore the evolution matrix for v{t) is E{t). On the other hand, the above definition can 
be written as 

Vi(t) = lq\t)Jiq(t), (54) 

where 

Jl = CTz, J2 = -CTx, J3 = 1, (55) 
and the cr's are Pauli matrices. Substituting Eq. (52) into Eq. (54), we have 

Vi{t)^y\0)[El{t)JiE,{t)]q{0). (56) 

Now any 2x2 matrix can be expanded in terms of the above Jj's and Jq = iay. Note that 
Jj are symmetric while Jq is antisymmetric, and the matrices E^{t)JiEq{t) are symmetric, we 
have 

El{t)JiE,{t)^aij{t)Jj. (57) 

It is easy to show that tT{JiJj) = 25^, and this yields aij{t) = ^tT[E^(t) JiEq(t) Jj]. Substituting 
into Eq. (56), we obtain Vi{t) — aij{t)vj{0). Therefore, Eij{t) — cnjit), that is 

E,,{t) = \ii[E\{t)J,E,{t)Ji]. (58) 

In other words, if the classical equation of motion is solved, which gives Eq{t), then E{t) can 
be obtained by simple algebraic calculations. This indicates a relation between our formalism 
and those of some previous authors, who find the time evolution operator of the Schrodinger 
equation by solving the classical equation of motion [6, 16]. However, our formalism, where e{t) 
plays the central role, is more convenient for the discussions of cyclic solutions and geometric 
phases, for example, in obtaining the necessary and sufficient condition (37). In practical 
calculations, it is usually more convenient to solve Eq. (17) directly than using the above 
relation. However, this relation is convenient for some general discussions. For example, when 
Eq{t) = ±1, it is obvious that Eij(t) = 6ij, 01 E{t) = 1. This will be useful below. 

Now we go into the main subject of this section. On the premise of Eq. (39), the time 
evolution operator is given by Eq. (40), where Q;(r) depends on the direction of Cq. If it 
happens that a(T)/7r is a rational number, then more cyclic solutions are available. Of special 
interest are the cases where q;(t) = 2NTr and q;(t) = (2A'" + l)7r. These will be discussed 
separately. 



A. q;(t) = 2iV7r 



In this case U{t) becomes a c-number. In fact, comparing Eqs. (40) and (4), we find that 
U{t) = Q{—2a{T), Co). If Q;(r) takes the above value, then U{t) = Q{—4:N7i, Bq), which has 
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been shown to be a c-number in Sec. II. The vahic of this number can be obtained by applying 
U{t) to a specific state, say, the ground state oi K ■ e^, which gives the result 

U{t) = e'^'^. (59) 

Several consequences can be deduced in this case. 

(1) All solutions are cyclic in the time interval [0,r], including nonnormalizable states, 
though we are only interested in normalizable ones. 

(2) Let Eq = (sinh,^o COS0O) sinh^o sin 00) cosh,^o)i where and 0o are arbitrary. One can 
choose tp{0) = iPq" = <5('Co) 0o)'^o as an initial state such that u{0) = /leo/2. In the state ip{t) 
with the above initial condition ip{0), we have u{t) = he{t)/2, where e{t) is the solution to 
Eq. (17a) with the initial condition eo, because u{t) and e{t) satisfy the same equation. Now 

u{t) = (V^(r), Ki;{T)) = (V^(0), t/t(r)Kt/(r)V^(0)) = (^^(0), K^) = "(0)- (60) 

Therefore 

e(T) = eo. (61) 

Because ^o and 0o are arbitrary, this means that the evolution matrix E{t) is a unit matrix at 
the time r: 

E{r) = 1. (62) 

In this case, we have obviously o"(r) = 1, that is, all three eigenvalues of E{t) are 1. The 
inverse is not true, however. In some cases we have three eigenvalues all equal to 1, but E{t) 
cannot be diagonalized, and is of course not a unit matrix (see Sec. VI). 

(3) Now we take eo, different from eo, as the initial condition for Eq. (17a), and use e{t) 
instead of e{t) in Eq. (34b), we obtain 

U(t) = exp[ih-^a{T)K ■ e^] , (63) 

where Q;(r) is given by Eq. (41), with e(t), ^(t) and 0(t) replaced by e(t), ^{t) and (f){t), 
repectively. This must be equal to that in Eq. (59), however. Thus we should have Q;(r) = 
27V7r, and N — N must be an even integer. Actually, we will show that q;(t) = q;(t), or N — N. 

Consider two initial unit vectors eo and eo, whose difference Sbq = Cq — Cq is infinitesimal 
(then deo-BQ = 0). The difference in Q;(r), according to Eq. (41), must be infinitesimal because 
the difference in e(t), and thus C,{t) and 0(t) are all infinitesimal. Therefore a(r) and thus 
Q;(r)/7r are continuous functions of Cq. Now that a{T)/7T can take only on even integers, an 
obvious consequence is that Q;(r) is a constant, independent of ^o and ^q. 

(4) Consider a cyclic solution in [0, r] with an arbitrary initial condition ■0(0) which is 
normalizable. The average value of K in ^(0) is denoted by w(0) as before. Since u^(0) > 
ft^/A > 0, we define Uo — ■\Ju'^{Q)/h which is dimensionless, and introduce 

So ^u{Qi)/huo. (64) 

Obviously, Cq = 1, and 693 > because ^3(0) > 0, thus this Cq can be used as the initial 
condition in Eq. (17). At later times, e{t) = u{t)/huo. The dynamic phase is 

p = -h-^ I uj{t)u{t) ■ n^(i) dt. (65) 
Jo 

Though Q;(r) = 2N7i is independent of Co, we must take the one given by Eq. (64) such that 
the second term in Eq. (41) can be related to the dynamic phase above. Then 



(3 = uo[a{T) + A^g] = uo{2Ntt + A^J, 



(66) 



13 



where A^g is calculated by substituting the above e(t) into Eq. (44). Because of Eq. (59), the 
total phase change is 5 = Nn, mod 27r. Therefore the geometric phase ^ — 5 — (3 turns out to 
be 

7 = -MoA^g - {uo - |)2iV7r, mod 2tt. (67) 

Here the first term is the familiar one, but an extra term appears, which depends on the initial 
condition. It vanishes (mod 27r of course) when mq — 1/2 is an integer, especially when the 
initial state is an eigenstate oi K ■ (it cannot be an eigenstate oi K ■ with some other Cq 
since otherwise u(0) would point in the direction of Cq) such that u(0) = knheo and Uq — kn- 
In the latter case it reduces to Eq. (43) as expected. 

(5) We have seen that Eq. (39) plus the condition Q;(r) = 2A^7r leads to Eq. (59), and as a 
result, all solutions arc cyclic in the time interval [0,r]. If e(t) is complicated, however, it is 
not convenient to use the above criterion because it may be difficult to calculate Q;(r). Thus 
some other convenient criterion is of interest. Now we give a necessary and sufficient condition 
for Eq. (59): 

E,(t)^1 ^ U(r)^e'''\ (68) 

First, suppose that U{t) — For arbitrarily given values Xq and po, it is easy to find 

an initial state ip{0) such that a;(0) = Xq and p(0) = Pq. It is obvious that x{r) = x{0), 
p(t) = p{0). Since x{0) and p{0) are arbitrary, we obtain Eq{r) = 1. Second, suppose that 
Eg{T) = 1. Then we have E{r) = 1 as mentioned below Eq. (58), and e(r) = Bq for any Bq. 
Now we choose Bq = (0,0, 1), and have from Eq. (40) U{t) = exp[i/i~^a(r)i^'3]. On the other 
hand, Eg{T) = 1 leads to x{t) — x{0) and p(r) = p(0) for an arbitrary ■ip{0), that is 

{m,u\T)xu{T)m) = {m,^m), {m,u\T)pu{T)m) = {m,pm)- m 

Since -0(0) is arbitrary, we should have U^{t)xU{t) — x and U^{t)pU{t) = p. From Eq. (14), 
this is valid only when q;(t) = 2A^7r, which leads to U{t) — d^'^. 



B. a(T) = {2N + l)7r 

In this case, U{t) = ex.p[m-^a{T)K ■ e^] = Q{Co,(l)o)exp[m-\2N + l)TrK3\Q''{Co,(l>o), 
where we have used Eq. (11) in obtaining the second equality. As pointed out in Sec. II, 
exp[ifi~^(2A'" + l)7rK3] commutes with K, and thus commutes with Q{$,o, (po), we have 

U{t) = exp[m-\2N + l)7rK ■ eg] = exp[in-^(2iV + 1)71X3] = e'^'' expiih-^K^). (70) 

Several consequences similar to those in subsection VA can be deduced in this case. 

(1) All normalizable solutions with definite parity arc cyclic in the time interval [0,r]. In 
fact, an initial state with definite parity can be expanded as 

00 00 

^"^(0) = X^'^2nV'2n, ^" (0) = ^ a2n+l^/'2n+l , (71) 
n=0 n=0 

where the superscript + (— ) indicates even (odd) parity, and time evolution does not change 
the parity of an initial state since U (t) only involves K and K is quadratic in x and p. It is 
easy to see that 

V'^(r)=exp(i5±)V'^(0), (72) 
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where 

5± = ±(A^ + i)7r, mod 27r. (73) 

(2) Repeat the discussions of the second point in subsection VA. Though f/(r) is not a 
c-number now, it commutes with K. Thus Eq. (60) is stiU vahd, and so is Eq. (62). 

(3) As before, we take Bq, different from Bq, as the initial condition for Eq. (17a), and 
obtain 

U{t) = exp[in-^a(T)ii: • eg] = Q(|o, 4>o) exp[ih-'a{T)Ks]Q\io, 4>o), (74) 

This must be equal to that in Eq. (70), however. Because exp[ih^^{2N + l)7ri^3] commutes 
with Q{Co:4>o): we have exp[ih~^a{T)K2] = exp[i^~^(2iV + l)7ri^'3]. This means that Q;(r) = 
(27V + l)7r, and N — N must be an even integer. By arguments similar to those in subsection 
VA, wc can conclude that q;(t) = a(r). 

(4) Consider a cyclic solution in [0, r] with a normalizable initial state ip{0) which is of 
definite parity. We define Bq as in Eq. (64) and use it as the initial condition in Eq. (17), then 
e(t) = u{t)/huQ. The dynamic phase is of the form in Eq. (65). As before, we use the above 
B{t) to calculate q;(t), then the second term in q;(t) can be related to the dynamic phase, and 

(5 = uo[a{T) + A^g] = uo[{2N + l)7r + A^g], (75) 

where A^g is calculated by substituting the above e(i) into Eq. (44). The total phase change 
has been given in Eq. (73). Therefore the geometric phase turns out to be 

7± = -MoA^g - {uq t i)(2A^ + l)7r, mod 27r. (76) 

Here an extra term appears once again, which depends on the initial condition. It vanishes 
(mod 27r of course) when mq =F 1/2 happens to be an even integer. For example, if the initial 
state is an eigenstate oi K ■ Bq (then it is of definite parity) such that tt(0) = knhBo and 
uq = km then the extra term vanishes as expected. 

(5) As in subsection VA, we give a necessary and sufficient condition for Eq. (70) which 
may be more convenient: 

E^{t) = -1 ^ U{t) = exp[ih-\2N + 1)71X3]. (77) 

First, suppose that U{t) — exp[i^~^(2A^ + l)7rKs]. Using Eq. (14), it is easy to show that 

U\t)xU{t) = -X, U\t)pU{t) = -p. (78) 

For arbitrarily given values Xq and po, it is easy to find an initial state ip{0) such that x(0) = Xq 
and p(0) = pq. The above equation yields x{t) = —x{0), p{r) = —p{0). Since x{0) and p(0) 
are arbitrary, we obtain Eg^r) = —1. Second, suppose that £'g(r) = —1. Then we have 
E{t) — 1, and b{t) — Bq for any Bq. Now we choose Bq — (0,0,1), and have from Eq. 
(40) U{t) = exp[i/i'~-^Q;(T)ir3]. On the other hand, -E'g(r) = —1 leads to x{t) = — a;(0) and 
p{t) = —p{0) for an arbitrary 'ip{0). Thus Eq. (78) must be true, and from Eq. (14), we have 
a{r) = {2N + l)7r, which leads to U{t) = exp[i;i-^(2iV + 1)ttKs]. 



C. Other cases 

Finally we briefly discuss the case where a{r)/7i is a rational number other than an integer. 
More specifically, let q;(t) = TqTi/sq where both tq and Sq are integers and prime to each other. 
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and So > 2. In addition to the denumerable set of cyclic solutions discussed before, we have 
also other ones in this case. For example, the initial conditions 

oo oo oo 

s=0 s=0 s=0 

gives cyclic solutions in the time interval [0,r]. In fact, it is easy to show that 'iP{t) = c"^-?/'(0) 
with S = Q;(r)/2 = ro7r/2so- In the above initial states, we have u{0) = hu^eQ where Uq = 
2sqs + 1/2 and s = ■^ksP- The dynamic phase is /? = mo[q!(t) +A9g]. Thus the geometric 
phase is 

7 = -UoA0g - 2roS7r, mod 27r. (80) 

We find once again that an extra term appears in the geometric phase. This term disappears 
when s happens to be an integer, especially when the above initial state involves only one 
term. 



VI. SOME EXAMPLES WITH EXPLICIT RESULTS 

In this section we calculate some simple examples where explicit analytic results are avail- 
able. The Hamiltonian involved in these examples may be elliptic, hyperbolic or critical. The 
motion of the particle exhibits various possible patterns. A normalizablc state can be regarded 
as a wave packet in the configuration space. The center of the wave packet is characterized 
by X and its velocity characterized by p, and their evolution is governed by the matrix Eg{t). 
If Eg{t) is finite at all times (for example, its elements are trigonometric functions of t), then 
the motion is said to be finite, because both x and p will remain finite at all times. Otherwise 
it is said to be infinite. The latter case contains still different patterns, for example, Eq(t) 
may increase with time or may be oscillating with increasing amplitude. Other quantities that 
characterize a wave packet are mainly the variances Ax = ^J {{x — xY) and Ap = \/ ((p — PY)- 
They can be regarded as the width of the wave packet in the coordinate and momentum space, 
respectively. Prom Eq. (58) we see that if Eq{t) is finite (infinite), then E{t) is essentially 
finite (infinite) as well. We will see that the time evolution of (Ax)^ and (Ap)^ are essentially 
governed by E{t), thus the width and the position of a wave packet exhibit similar patterns 
of motion in this system. In other words, if the wave packet moves in a confined region, then 
its width also remains finite, and the contrary is also true. In fact, we can define a new vector 
w{t) by 

w{t)^u{t)-v{t). (81) 

Since both u{t) and v{t) satisfy the same equation as e(t), it is obvious that w{t) also satisfy 
the same equation, and thus its time evolution is governed by the matrix E{t). On the other 
hand, it is easy to show that 

w, = \[{Axf - {Apfi = \[{Axf + {Ap)% (82a) 
W2 = -\{{x-x){p-p) + {p-p){x-x)). (82b) 

Therefore, 

{Axf = W3 + wi, {Apf = - wi, (83) 
and their time evolution is essentially governed by E{t) . 



16 



We will see that the pattern of motion does not have a simple relation with the nature 
of the Hamiltonian. For a Hamiltonian with a definite nature, say, elliptic, the wave packet 
may exhibit different patterns of motion if some parameter in the Hamiltonian takes different 
values. And, when the parameter goes through some critical value, the motion changes from 
one pattern to another. This is somewhat like phase transition, and seems not to be noticed 
before. 



A. n = constant 



The first case is the simplest one where n is a constant vector and u!{t) is an arbitrary 
function of time. Let 



(p{t) = [ u{t') dt, 
Jo 



(84) 



and if e{t) = e[ip{t)] obeys 

e'(^) = -2n^ X e^i^) (85) 

where the prime indicates derivative with respect to ip, then Eq. (17a) is satisfied. This is 
similar to Eq. (5) and can be solved in the same way. The result depends on the nature of n, 
and will be given separately. 

(1) If = 1, the solution reads 

e{t) = [bq — (eo • n^)n] cos 2cp + Bq x -nP sin + (cq • n^)n. (86) 

To obtain the evolution matrix E{t), we should write down the components of n: 

n = (sinh ^„ cos 2(^„, sinh ^„ sin 2<^„, cosh ^„) , (87) 

where ^„ and are arbitrary. The matrix E{t) is somewhat complicated, but it can be written 
as the product of several simple matrices 



E{t) = R{iPn)S{U, l)Wi'\if )S-\U, m~\^n), 



where 



^cos 2(p — sin 2ip 0^ 



sin 2(p cos 2(p 


and the matrices R and S, which will be repeatedly used below, are defined by 



V 



(88) 



(89) 



^cos2ipn — sin2(/7„ 0^ 



Ri'Pn) 



sin 2(pn cos 2ipn 
1 



^ecosh^n sinh^„ 



1 

sinh^„ ecosh^j 



(90) 
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where e = ±1. The eigenvalues of E{t) can be easily found to be {1, e^^''', e"'^*^}. The eigen- 
vector corresponding to the eigenvalue 1 is ry = n. This can be easily seen from Eq. (86), 
where Bq — n leads to e(i) = n at any time t. Because — 1, there exist normalizable cyclic 
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solutions at any time interval [0, r]. However, most of these cyclic solutions are trivial ones. 
Actually, the time evolution operator can be obtained by direct integration in this case: 

U{t) = exp(-in-V-'^ • ^^)- (91) 

If ip{0) is an eigenstate of K ■ n^, then ilj{t) is different from ip{0) only by a phase factor, thus 
it is cyclic in any time interval but is trivial. Only when ^p{t)/t: takes rational numbers that 
we have nontrivial cyclic solutions in [0,r]. 

The evolution matrix Eq{t) can be similarly found to be 

E,{t) = R,{^n)w!^^{^)R-\^r.), (92) 

where 

(n-\ I COS </? exp(£„) sin (/? , 

\ — exp ( — ^„) sin (/? cos </? 



and 



cos<^„ smipn \ 

Rqi^n) =1 ■ (94) 

sin cos (/9„ / 

Obviously, when iy9(r) = —Nrc, Eqij) — (—1)"^ and E{t) = 1. This includes the two most 
interesting cases discussed in Sec. V. 

The time evolution operator (91) can also be obtained by our method. Taking the simplest 
solution e{t) — n, we have ^{t) = = Cn, 0(0 — (l>o — 2</7n, and a{t) — —^p{t). When these 
are substituted into Eq. (34b), we obtain the above result. Moreover, the above condition 
ip{t) = —Nn for -E'q(r) = (— 1)*^ is equivalent to q;(t) = Ntt, as expected. 

To verify Eq. (67) by explicit calculations, we take n = (0,0,1), then U{t) = 
exp[—ih~^ip{t)Ks]. When (^(r) = —2Ntt, we have U{t) = e'^'^. In the time interval [0,r], 
all solutions are cychc with 6 — Nn, mod 27r. Let the initial state have w(0) = huoCo, 
where Cq = (sinh cos 0o,sinh^o sin cosh ^o)- Then u{t) — huoe{t), where e{t) — 
(sinh ^0 cos(2(/7+0o) , sinh sin(2</7+0o) , cosh ^o) ■ It is then easy to find that /3 = 2Nt^uq cosh ^o, 
and 7 = A'^7r(l — 2Mocosh,^o)- On the other hand, the Hannay angle can be found to be 
A6'g = 2A^7r(cosh^o ~ !)• With these results it is easy to see that Eq. (67) is valid. 

(2) If = —1, the solution reads 

e{t) = [eo + (eo • n^)n\ cosh 2(/7 + x sinh 2(/7 - (cq • n^)n. (95) 
To obtain the evolution matrix E{t), we should write down the components of n: 

n — (cosh cos 2(^„, cosh ^„ sin 2(^„, sinh ^„) , (96) 
where ^„ and are arbitrary. The matrix E{t) is given by 

E{t) = R{ipn)S{U l)W^^\ip)S-\U m~\^n), (97) 

where 



^10 ^ 

cosh 2ip — sinh 2ip 
yO — sinh 2ip cosh 2ip j 



(98) 
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The eigenvalues of E{t) can be easily found to be {1, e^*^, e"^*^}. The eigenvector corresponding 
to the eigenvalue 1 is r/ = n. Because = —1, there exists no normalizable cyclic solution 
at any time interval [0,r]. The evolution matrix Eq{t) can be found to be 



(99) 



where 



(100) 



cosh (f exp (^„) sinh cp 
^exp ( — ^n) sinh if cosh (p 

In this case Eq(r) cannot take the value ±1, otherwise one should obtain Wq^(ip) — ±1. The 
latter is obviously impossible, unless </?(t) returns to but this is not of interest. 
(3) If = 0, the solution reads 



e{t) = 2ip'^{eo ■ n^)n + 2(pel x + eo- 



:ioi) 



Up to a constant factor, which may be absorbed in uj{t), the components of n may be written 
as: 

n = (cos 2ipn, sin 2ipn, 1), (102) 
where </?„ is arbitrary. The matrix E{t) is given by 

E{t)^R{iPn)wi;'\ip)R-\iPn), 

where 



(103) 



/l _ 2(^2 _2(^ 2(^2 \ 
2(/7 1 -2(/? 

,2 r,,„ -1 I r,,J2, 



(104) 



The eigenvalues of ii^(t) can be easily found to be {1, 1, 1}, but there is only one eigenvector 
T] — n a (p 0. Thus we find an example where all eigenvalues are 1 but E{t) cannot be 
diagonalized, let alone be a unit matrix. Because = 0, there exists no normalizable cyclic 
solution at any time interval [0, r]. The evolution matrix Eq{t) can be found to be 



E,{t) = R,i^r.)W^i^)Rq\^^), 



where 




(105) 



(106) 



In this case -Eq(r) cannot take the value ±1 either, unless </7(t) returns to but this is not of 
interest. 

From the above results we see that the motion is finite when the Hamiltonian is elliptic and 
infinite (assume that ip{t) can reach arbitrarily large values, for example, Lp{t) is proportional 
to t) in the other cases. But note that E{t) and Eq{t) arc polynomials of p in the critical case 
and are exponential functions of it in the hyperbolic case. If n is time dependent, however, the 
situation is not so simple. We will see that for a Hamiltonian of a definite nature, all patterns 
of motion are possible. 
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B. n = (ni cos 2(^, ni sin2(^, na) 

Here rii and are constants, satisfying = n1 — nf = ±1 or 0, and (/? = Lp{t) is an 
arbitrary function of t with ^^(O) = 0. For convenience, we assume henceforward that ip{t) 
increases with t monotonically and goes to infinity when t does. If u!{t) — X(p{t) where A is a 
constant, analytic solutions are available. 

By making a time-dependent unitary transformation U{t) — exp[—ih~^(fiK3HJ{t), it is not 
difficult to obtain 

U{t) = exp(-i^-Vi^3) exp{-ih-^A(pK ■ n^), (107) 
where A > and nj = ±1 or 0, and 

An/ = (Ani,0, Ana - 1). (108) 

If rij = 1, wc sec that cyclic solutions are available in [0,r] if the initial state is an eigenstate 
oi K -n^jr and v?(r) = Nn, and more general cyclic solutions are available in some time interval 
if A happens to be a rational number. However, for other time interval, or when nj — —1 or 0, 
the above time evolution operator, though explicit, docs not tell much about cyclic solutions. 
Consequently, it is necessary to study the matrix E{t) and Eq{t). 
If e{t) = e[cp{t)] obeys 

e'{^) = -2Xn%^) X e^{ip) (109) 

then Eq. (17a) is satisfied. Now we make a time-dependent linear transformation, written in 
column vector form: 

e{ip) = Riip)f{ip), (110) 

where the matrix R is defined in Eq. (90), but now the independent variable (f is time 
dependent. The reduced equation for /(</?), written in ordinary vector form, reads 

f{ip) = -2An] X n^). (Ill) 

Now n,/ is a constant vector, so the reduced equation is easy to solve. The form of the solution 
depends on the nature of n/, and several cases should be treated separately. 

(1) If = 1 and A G (—00,723 — |ni|) U {n^ + |ni|,+cxo), or -n? = —1 and A G (— ns — 
~f^3 + I^iDj or = and Ana < 1/2, then (Ans — 1)^ > (Ani)^. In this case we define 

A = y/iXns - 1)2 - (Ani)2, cosh e„ = | Ana - 1|/A, sinh^n = Am/A, (112) 

then 

n/ = (sinh^„,0,ecosh^„), (113) 

where e = e(Ana — 1) is the sign function of Ana ~ 1; that is, it equals 1 (—1) if Ana — 1 > 
(< 0). The solution is 

E{t) = Riip)SiU,e)wi'\^)S-\U,e), E,{t) = R,iip)W^'^iip), (114) 

where 



( cos Ikip - sin 2A99 0^ 
sin 2A99 cos 2A(y9 

V " ° V 



COS A(/7 e exp(e^„) sin Ac/? 

- e exp ( — ) sin Ac/? cos A(/7 



(115) 
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and the matrix Rq is defined in Eq. (94), bnt now the independent variable (p is time dependent. 
In this case the motion is finite, no matter the Hamiltonian is eUiptic, hyperbohc or critical, 
as long as A belongs to the appropriate region. 

The eigenvector of E{t) corresponding to the eigenvalue 1 is (not normalized) 

ri{t) — (esinh^„ sin A</7 cos </?, esinh^„ sin A(/7 sin </?, cosh ^„ sin A</7 cos </? + e cos A(/7 sin (^). (116) 

If sin(/7 = but sinA0 7^ at the time r, or sin Ac/? = but sincf) ^ Q (if both are zero, then 
E{t) — 1), then rf{T) > and there are cyclic solutions in the time interval [0, r]. In the 
general case, we have 

'if'{t) = (sinh^ ^„ sin^ A(^ + 1) sin^[(/9 + earctan(cosh,^„ tan A(^)] — sinh^ sin^ A<y9. (117) 

This may be either positive or negative. If at the time r the quantity in the square bracket is 
close to (A^ + l/2)7r, then 'r7^(r) > 0, and there are cyclic solutions in the time interval [0,r]. 
If the quantity is close to Nn, and sinA^? is not too small, then r;^(r) < 0, and there is no 
cyclic solution in the time interval [0, r]. 

If A is a rational number, then there exist some r such that Eq{T) — ±1 and E{t) — 1. 
This is the case where more general cyclic solutions are available. 

The time evolution operator can be easily obtained by our formalism in this case. We take 
the simple solution e{t) = (e sinh^„ cos 293, esinh^„ sin 293, cosh ^„) and using Eq. (34a), after 
some operator algebras, we arrive at the result (107). For the subsequent cases, however, our 
formalism is not convenient in calculating the time evolution operator, since it is not easy to 
find a solution that is sufficiently simple and satisfies Eq. (17c). Instead, the time-dependent 
unitary transformation that leads to the result (107) is convenient in this respect. However, 
as pointed out before, the result (107) does not tell much about cyclic solutions. 

(2) If = 1 and A G (na — na + |ni|), or = —1 and A e (—00, —n^ — \ni\) U (— na + 
|ni|, +00), or = and Xn^ > 1/2, then (Ana — 1)^ < (Ani)^. In this case we define 

A = ^/(ArHF-^TAn^^^, cosh^„ = |Ani|/A, sinh^„ = (Arii - 1)/A, (118) 

then 

Uf = (ecosh^„,0,sinh^„), (119) 

where e — e{Xni) . The solution is 

E{t) = R{^)S{U,e)W^'\^)S-\Ue), Eq{t) = Rq{^)W^1{^), (120) 



where 



^10 ^ 

cosh 2A.LP — sinh 2A.ip 
^0 — sinh 2A(/7 cosh 2A</7 ^ 



cosh A(/7 e exp (e^n) sinh A</7 

e exp ( — ) sinh A93 cosh Ay? 



(121) 

In this case the motion is oscillating with exponentially increasing amplitude, no matter the 
Hamiltonian is elliptic, hyperbolic or critical, as long as A belongs to the appropriate region. 
The eigenvector of E{t) corresponding to the eigenvalue 1 is (not normahzed) 

ri{t) — (e cosh ^„ sinh A.ip cos </?, e cosh ^„ sinh A.ip sin (/?, sinh ^„ sinh Ac/? cos </? + cosh A(/j sin Lp) . 

(122) 
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If sin ip = but (/? 7^ at the time r, then t}^{t) < and there is no cychc solution in the time 
interval [0, r] . In the general case, we have 

r]^{t) — (cosh^ sinh^ Acp + 1) sin^[</7 + arctan(sinh ^„ tanh Acp)] — cosh^ ^„ sinh^ Acp. (123) 

This may be cither positive or negative. If at the time r the quantity in the square bracket is 
close to (A^ + l/2)7r, then r7^(r) > 0, and there are cyclic solutions in the time interval [0, r]. 
If the quantity is close to Ntt, and ip is not too small, then ?7^(r) < 0, and there is no cyclic 
solution in the time interval [0, r] . In the present case there exists no r such that (/?(t) ^ 
and Eq{T) = ±1. 

(3) If = 1 and A = ns ± \ni\, or = —1 and A = —n^ ± or = and Xn^ — 1/2, 
then (Ans — 1)^ = (Ani)^. Let e = (Ana ~ 1)/Ani = ±1. The solution is 



E{t) = R{^)D^'\^), E,{t) = RMW^Si^), (124) 



where 



( 1 - 2(Ani) -2eAni(/? 2e(Ani) ^ 

2e\n\kp 1 —2Xn\ip 

^-2e(Ani)V^ -2Ani(^ 1 + 2(Ani) 



^ 1 (1 + e)\ni^ 

y(l — e)Ani</7 1 



(125) 

In this case the motion is oscillating with polynomially increasing amplitude, no matter the 
Hamiltonian is elliptic, hyperbolic or critical, as long as A takes the appropriate value. 
The eigenvector of E{t) corresponding to the eigenvalue 1 is (not normalized) 

ri{t) — (Ani(/p cos </?, Ani</? sin (/?, eAni(/p cos </p + sine/?). (126) 

If simp = at the time r, then rf{T) = and there is no cyclic solution in the time interval 
[0, r]. In the general case, we have 

rj'^(t) = (A^ni<^^ + 1) sin^[(/7 + e arctan(Ani(/7)] - A^n^c^^. (127) 

This may be either positive or negative. If at the time r the quantity in the square bracket is 
close to {N + l/2)7r, then T7^(t) > 0, and there are cychc solutions in the time interval [0, r]. 
If the quantity is close to Nn, and ip is not too small, then T7^(t) < 0, and there is no cyclic 
solution in the time interval [0, r]. In the present case there exists no r such that (p{t) ^ 
and Eq{T) = ±1. 

It should be noted that all the three eigenvalues of E{t) are 1 when ip{t) = Nn, and all the 
two eigenvalues of Eq(t) are 1 when (p{t) = 2Ntt, but they are not unit matrices unless N — 0. 

Prom the above results we see that for a Hamiltonian of a definite nature, say, elliptic, the 
motion may have all possible patterns, depending on the value of A. When A < ns — |ni| or 
A > ns + |ni| the motion is finite. When ns — |ni| < A < ns + |ni|, the motion is oscillating 
with exponentially increasing amplitude. At the two critical values A = ns ± |ni| the motion is 
oscillating with polynomially increasing amplitude. Something like phase transition happens 
here. When A goes through the critical values, the motion changes from one pattern to another, 
and at the critical values the motion exhibits an independent pattern. Similar situation will 
be seen in the next subsection. 

A common feature of the three different cases above is that in any time interval (tq, +00) 
with To > one can always find some r such that cyclic solutions are available in [0, r]. 
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C. n = (ni, na sinh2(^, 77,3 cosh 2<^) 

As before, rii and are constants, satisfying = n| — = ±1 or 0, and = ip{t) 
has the same property as in subsection VI B. If ijj{t) = X(p{t) where A is a constant, analytic 
solutions are available. The time evolution operator can be obtained in a way similar to that 
in subsection VI B. The result is 

U{t) = exp(-i^-Vi^i) exp{-ih-^A(pK ■ nj), (128) 

where A > and rij — ±1 or 0, and 

An/ = (Ani + 1,0, Ans). (129) 

Compared with the case in subsection VI B, here still less about cyclic solutions can be seen 
from the above result. Consequently, it is necessary to study the matrix E(t) and Eg(t). 

If e{t) = e[ip{t)] obeys an equation of the form (109), then Eq. (17a) is satisfied. As before 
we make a time-dependent linear transformation, written in column vector form: 

e((^)=r((^)/((^), (130) 

where the matrix T, and another matrix Tg that will be used below, are defined as 

A \ 



cosh2(/7 sinh2(/7 
yO sinh2(/7 cosh2(/7^ 



^ ^ , coshv? -sinh(/9 , 

TM =1 I • (131) 

sinh (p cosh (p 



The reduced equation for /((/?), written in ordinary vector form, has the form of Eq. (HI), 
but where Anj is given by Eq. (129). As before, several cases are to be treated separately. 

(1) If = 1 and A e (— oo,ni — |n3|) U (ni + |n3|,+oo), or = —1 and A e (— ni — 
1^3 1; — ""-1 + 1^3 1); or Ti^ = and Xrii < —1/2, then (An3)^ > (Ani + 1)^. In this case we define 

A = V(An3)2 - (Ani + 1)2, coshen= |An3|/A, sinh^n = (Arii + 1)/A, (132) 

then Uf has the form in Eq. (113), but e = e(An3). The solution is 

E{t) ^T{ip)S{^n,e)wi'\ip)S-\^n,e), Eg{t) ^ T^W^'^ {^) , (133) 

where w!i!\(p) and Wg'^{(p) have been given in Eq. (115). In this case the motion is oscillating 
with exponentially increasing amplitude. 

The eigenvector of E{t) corresponding to the eigenvalue 1 is (not normahzed) 

ri{t) = (e(sinh ^„ cosh ip sin Aip — sinh ip cos Aip) , cosh ^„ sinh ip sin Aip, cosh ^„ cosh ip sin Aip) . 

(134) 

If smA(p — but (p at the time r, then r]'^{T) < and there are no cyclic solution in the 
time interval [0,t]. In the general case, we have 

77^(t) = cosh^ sin^ A9? — (cosh^ cosh^ ~ 1) sin^fA^? — arctan(tanh99/ sinh,^^)]. (135) 

This may be either positive or negative. For large (f the second term is large as long as the 
quantity in the square bracket is not close to Nt:, then r]'^{T) < 0, and there is no cyclic 
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solution. However, if (p{t) is close to lpq where (po is the root of tanA(/9 = tanhLp/ sinh^„, the 
second term is very small and 77^ (r) > 0, then we have cyclic solutions in the time interval [0, r]. 
The above transcendental equation has infinitely many roots, therefore in any time interval 
(tq, +00) with To > we can always find some r such that cyclic solutions are available in 
[0, r]. In the present case there exists no r such that 9?(r) 7^ and -Eg(r) = ±1. 

(2) If = 1 and A e (rii — Iriaj, ni + |n3|), or = —1 and A e (—00, —rii — \ns\) U (— rii + 
jrial, 00), or = and Xrii > —1/2, then (Ana)^ < (Arii + 1)^. In this case we define 

A = V(Ani + 1)2 - (Ana)^, cosh^^ = |Ani + 1|/A, sinh^^ = Ang/A, (136) 

then Uf has the form in Eq. (119), but e = e(Ani + 1). The solution is 

E{t) = T{^)S{U,e)W^_'\ip)S-\^n,e), E,{t) = T,{^)W^}{^), (137) 

where w'^\ip) and W^}{ip) have been given in Eq. (121). In this case the motion is exponen- 
tially infinite. 

The eigenvector of E{t) corresponding to the eigenvalue 1 is (not normalized) 

'q{t) — (e cosh ^„ cosh (p sinh K(p — sinh (/? cosh Ac/?, sinh ^„ sinh (/? sinh A(/?, sinh ^„ cosh (/? sinh A(/?) , 

(138) 

which gives 

rf{t) = sinh^ sinh^ A99 — (cosh ^„ cosh 9? sinh A99 — e sinh 9? cosh A9?) 2. (139) 

If e = —1, it is obvious that the second term is larger and rf{T) < 0, and thus there is no 
cyclic solution in any time interval [0, r]. If e = 1, the above result can be recast in the form 

rfit) ~ sinh^ ^„ sinh^ A9? — (sinh^ ^„ cosh^ + 1) sinh^[A93 — Artanh(tanh ip/ cosh^„)]. (140) 

In this case, cyclic solutions are possible. For example, when A < 1, the transcendental 
equation tanh(/9 = cosh ^„ tanh A(/9 has one root (pQ. When ^{t) is close to (/9o, the second term 
is very small and rfir) > 0. Then we can have cyclic solutions in [0, r]. Another possible case 
is </?(t) <^ 1, which also leads to 'q'^ir) > if exp(|^„|) > A. Then we also have cychc solutions 
in [0, r]. However, when (p{t) is large, it is easy to see from Eq. (139) that tj'^{t) < regardless 
of the parameters ^„ and A. Then there is no cyclic solution in [0, r]. In other words, no initial 
state can return to itself at a large time. This is rather different from the last case. As before, 
there exists no r such that (/^(r) 7^ and Eq{T) = ±1. 

(3) If = 1 and X — ni± \ns\, or n? — —1 and A = — ni ± \ns\, orn^ = and Ani = —1/2, 
then (Ans)^ = (Ani + 1)^. Let e = (Ani + 1)/An3 = ±1. The solution is 

E{t) = T{^)D^^\^), E,{t) = T,{^)W^^{^), (141) 

where D^'^\ip) and W^q{(p) can be obtained from D''^\ip) and W'^{ip) in Eq. (125), respec- 
tively, by substituting eXnz for Arii. In this case the motion is also exponentially infinite, but 
different from the last case. Indeed, both w'^\ip) and W^}{ip) in Eq. (137) involve exponential 

functions of </?, while D^'^^ip) and Wqo^(</') involve polynomials of it. 

The eigenvector of E{t) corresponding to the eigenvalue 1 is (not normahzed) 



rj{t) = {eXn^ip cosh if — sinh ip, Xn^ip sinh ip>, Xn^tp cosh 93), 



(142) 
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which gives 

ri'^{t) — {Xn^^pf' — (eAn3(/7COsh(^ — sinh(^)^. (143) 

If eAris < 0, it is obvious that the second term is larger and rfir) < 0, and thus there is 
no cychc solution in any time interval [0,r]. If eXn^ > 0, cyclic solutions are possible. For 
example, when |An3| < 1, the transcendental equation tanh(/9 = eXn^ip has one nonzero root 
i-Pq. When ^{t) is close to ipo, the second term is very small and ffij) > 0. Then we can have 
cychc solutions in [0, r]. Another possible case is (p{t) -C 1, which also leads to r7^(r) > and 
we have cyclic solutions in [0, r]. However, when (p{t) is large, it is easy to see that »7^(t) < 
regardless of the parameters and A. Then there is no cyclic solution in [0,r]. In other 
words, no initial state can return to itself at a large time. This is similar to the last case. As 
before, there exists no r such that (^(r) ^ and Egi^r) = ±1. 

In this subsection we also see that different values of A correspond to different patterns of 
motion, for a Hamiltonian of a definite nature. But the cases are rather different from those 
in subsection VI B. In some of the cases here, there is no cyclic solution in any time interval. 



D. n = (ni, n3 cosh2(^, n3 sinh2(^) 

Here rii and ns are constants, satisfying — —n\ — n\ — —1. This is different from the 
previous cases, since cannot take 1 and 0. = ip{t) has the same property as before. 
If uj{t) = X(f{t) where A is a constant, analytic solutions are available. The time evolution 
operator is given by Eq. (128) but now 

An/ = (Ani + 1, An3,0). (144) 

The equation for e{t) can be solved in a way similar to that in subsection VI C. We define 

A = y/{Xni + 1)2 + (An3)2, cos2(^„ = {Xrii + 1)/A, sin2(^„ = An3/A, (145) 

then Uf — (cos2(/7„,sin2(/?„,0). The solution is 

E{t) = T{ip)R{ipr,)w['\ip)R-\ipr,), E,{t) = TMR,{ipr,)W^''\ip)R-\ipr,), (146) 
where w!i\(p) is given in Eq. (121) and 

, ,s / cosh Acp sinh A(p \ 

Wt\^) = ^ (147) 

\ sinh A(p cosh A(p I 

Obviously, the motion is exponentially infinite. 

The eigenvector of E(t) corresponding to the eigenvalue 1 is (not normahzed) 

ri{t) = (cos 2ipn cosh ip sinh Ap — sinh ip cosh Ap, sin 2p>n cosh p> sinh Ap>, sin 2p>n sinh p> sinh Aip) , 

(148) 

which gives 

(t) = — sin^ 2(^„ sinh^ A(^ — (cos 2(^„ cosh sinh A(^ — sinh <^ cosh A(^)^. (149) 

It is obvious that t}'^{t) < if </?(t) 7^ 0, and there is no cyclic solution in any time interval 
[0, r]. In other words, in this case no initial state can return to itself at a later time. As before, 
there exists no r such that (^(r) 7^ and Eg^r) = ±1. 
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VII. SUMMARY 

In this paper wc develop a method for solving the Schrodinger equation of the generalized 
time-dependent harmonic oscillator. This method, though not always convenient in practical 
calculation of the time evolution operator, is very suitable for the study of cyclic solutions and 
geometric phases. We concentrate our attention on Hamiltonian of general time dependence 
and cyclic solutions in the time interval [0, r] with an arbitrarily given r. A necessary and 
sufficient condition for the existence of cyclic solutions in such time intervals is given. There 
may exist some time interval in which more solutions are cyclic. This includes several cases 
among which two are of more interest. In one of these cases all solutions are cyclic, and in an- 
other all solutions with definite parity are cyclic. Criterion for the appearance of such cases are 
given. The proportional relation between the nonadiabatic geometric phase and the classical 
Hannay angle is reestablished. However, this holds only for cyclic solutions with special initial 
conditions. For more general cyclic solutions encountered in the above cases, the nonadiabatic 
geometric phase contains in general an extra term in addition to the one proportional to the 
classical Hannay angle. Several examples are studied where the Hamiltonians are relatively 
simple and analytic solutions are available. In these examples the existence of cyclic solutions 
are discussed in detail. They exhibit many possibilities: (1) cyclic solutions are available for 
all r; (2) available for some r and such r may be found in any interval (tq, -|-cxo) with tq > 0; 
(3) available for some r but such r exists only in some finite interval (0,ro); (4) not available 
at all. From the point of view of wave packets, the motion in these examples also exhibits 
various patterns. For a Hamiltonian of a definite nature, say, elliptic, several different patterns 
of motion are possible, depending on the value of some parameter in the Hamiltonian. There 
exists some critical value of the parameter at which some kind of phase transition happens. 
When the parameter goes through it, the motion changes from one pattern to another, and at 
the critical value itself the motion has an independent pattern. 
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APPENDIX 

Our approach to the solution of the Schrodinger equation with time-dependent Hamiltonians 
was first used in our previous paper [36] to the system of general spin moving in an arbitrarily 
varying magnetic fields, and is further developed in this paper. The Hamiltonian for the spin in 
a magnetic field is an element of the SO (3) algebra while that for the time- dependent harmonic 
oscillator is one of the SO (2,1) algebra. In this appendix we briefly discuss how to extend the 
formalism to a system where the Hamiltonian is an element of a more general Lie algebra. 

Consider a Lie algebra with the generators Ka (a = 1, 2, . . . , /), satisfying the commutation 
relation 

[Ka,K,]^iWK,, (A.l) 
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where fab^ are the structure constants. Define 

5a(, = /a//6/ = tr(/„A), (A.2) 

where /„ is a matrix whose element is defined by {fa)b^ — fab^- It is obviously symmetric: 
Qba = gab- FoT a semisimple Lie algebra, the matrix gab is not singular, and thus invertible. 
Its inverse matrix is denoted by g"''^, satisfying g°'^gd} = ^"'b- The matrices gab and g""^ acts 
like metric tensors and can be used to lowering and raising vector indices. For the SO (2,1) 
algebra with the generators given in Eq. (2), one should find gab — 8diag (1, 1, — 1) and 
g"^ — (1/8) diag (1, 1, —1). Thus the metric used in the main text is different from this by a 
constant factor. 

For simplicity we consider a semisimple Lie algebra of rank 2, say, the SU(3) algebra. It 
contains two mutually commuting operators which generates the Cartan subalgebra. Thus 
there are two Casimir operators, given by 

C2 = g'^'KaKb, Cs = h'^'^KaKbR,, (A.3) 

where h"'^'^ is obtained by raising the indices of habc, and the latter is defined by 

habc^tr{fafbfc)- (A.4) 

Now consider a physical system whose Hamiltonian is an element of the above semisimple 
Lie algebra. More specifically, 

H(t) = nu;%t)Ka, (A.5) 

where u!"'{t) are time-dependent frequency parameters, and the generators Ka are dimension- 
less. We define a Z-component vector u{t) — {ui{t) , U2{t) , . . . ,ui{t)) by 

Ua{t) = {Ka) = {m,Kam), (A.6) 

where ip{t) is an arbitrary state of the system, that is, a solution to the Schrodinger equation. 
It is easy to show that it satisfies the equation 

Ua{t) = fab'oo\t)u,{t). (A.7) 

This is a system of / linear differential equations. Now we define a /-component vector e{t) — 
(ei(i), e2(i), . . . , e.i{t)) by the same equation, that is 

ea{t)^fab''uj\t)e,{t), (A.8) 

and a nontrivial (nonzero) initial condition. It is then straightforward to show that the two 
operators linear in Ka, defined by 

Li{t) = g''''ea(t)Kb, L^H) = h''''ea(t)eb(t)K, (A.9) 

are invariant operators, and they commute with each other. Therefore, the invariant operators 
can be obtained by solving a linear differential equation for e{t). It can be similarly shown 
that g'^ea{t)eb{t) and h'^'^ea{t)eb{t)ec{t) are also time-independent quantities. Thus among 
the I components of e{t) only / — 2 arc independent variables. We can parameterize e(t) 
in some way similar to that in Eq. (15). The independent parameters will be denoted by 
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The subsequent steps in obtaining the time evolution operator depend on the details of the 
Lie algebra. We can give only a rough sketch here. 

The crucial step is to find a unitary operator Q{t) = Q{i{t)) such that 

L,it) = Q{t)H^Q\t), L^it) = Qit)H2QHt), (A.IO) 

where Hi and H2 are generators of the Cartan subalgebra. It might be difficult to prove the 
existence of Q{t) in the general case. For a specific Lie algebra, however, one can try to find 
it by practical calculations like those in Sec. 11. Since Hi and H2 commutes with each other, 
they have a complete set of common eigenstates. These will be denoted by {ipnm}, satisfying 

Hll/Jnm = ^n'^nm, -f^2V'nm = l^m'^nm- (A-ll) 

Obviously, A„ and /x^ are also the eigenvalues of Li{t) and L2{t), respectively. 
Next we will show that if the initial state -0(0) of the system satisfies 

Li(0)V^(0) = Xnilj{Q), ^2(0)^(0) = iim^iQ), (A.12) 

then the state tlj{t) at later times will satisfy 

Li(t)^lj(t) = Xn^l^it), L2{t)i:{t) = IXmi^it). (A.13) 

We use here a different proof from the induction one employed in Sec. 111. It is easy to show 
that = [Liit) - \n]il^{t) satisfies 

ihdt^l^^^\t) = H{t)ij^^\t), (A.14) 

where the equation ihdtLi{t) + [Li{t)^H{t)\ — has been used. Consequently, il)'^^\t) — 
U(t)il}^^\Q) where U{t) is the time evolution operator of the Schrodinger equation. Because 
Eq. (A.12) leads to ip^^\Q) = 0, we have ip^^\t) = 0, and thus the first equation in (A.13) is 
obtained. The second one can be proven in a similar way. 
On account of the above results, we have 

i>{t) = exp[ianni{t)]Q{t)i>nm: (A. 15) 

where Q;„^(i) is a phase that cannot be determined by the eigenvalue equation. However, it 
is not arbitrary. By the requirement that ■?/'(t) satisfies the Schrodinger equation, it can be 
determined in terms of ^(t). Finally, one should manage to replace A„ appearing in Q;„rn(^) 
and anrni^) by Hi, and /^.^ by H2, and obtain 

jP{t) = Q{t) eMHHi, H2, i)]Q^(0)V(0). (A.16) 

Because all the operators in the above equation are independent of n and m, and because the 
set {ipnm} is complete, we obtain the time evolution operator 

U{t) = Q{t) exp[ia(i/i, H2, t)]Q^(0). (A.17) 

The time dependence in both Q{t) and a{Hi, H2.t) comes from ^{t) (and uj"'{t) of course). 
Therefore, the time evolution operator is also obtained by solving the linear different equation 
for e{t). 
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